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SUMMARY 


The purpose of this study was to evaluate the use of a computer -generated perspec- 
tive tunnel display for following a strongly curved flight path. The display was evaluated 
by monitoring pilot performance in a fixed -base simulator with the vehicle dynamics of a 
CH-47 tandem-rotor helicopter, linearized about a 33.5-m/sec (110 ft/sec) forward 
velocity and a level flight trim condition. Superposition of the predicted future vehicle 
position on the tunnel image was also investigated to determine whether, and to what 
extent, it contributes to better system performance (the best predicted future vehicle 
position was sought). 

Three types of simulator experiments were conducted: following a desired tra- 

jectory in the presence of disturbances; entering the trajectory from a random position, 
outside the trajectory; and detecting and correcting failures in automatic flight. Two 
other display configurations were considered for comparison: a tracking display which 

incorporates future trajectory information by displaying the predicted future error, and 
combined vertical and horizontal situation displays containing aircraft attitude informa- 
tion and trajectory map information, respectively. 

The tunnel display with superimposed predictor /director symbols was shown to be 
a very successful combination, which outperformed the other two displays in all three 
experiments. A prediction time of 4 to 7 sec was found to optimize trajectory tracking 
for the given vehicle dynamics and flight condition. Pilot acceptance of the tunnel plus 
predictor /director display was found to be favorable, and the time the pilot needed for 
familiarization with the display was found to be relatively short. 

The simulator study was conducted on the Langley Real-Time Simulation System 
using a fixed-base simulator cabin with helicopter controls. An Adage Graphics 
Terminal, the AGT/130, was programmed to generate the various display configurations. 
The Control Data CYBER 175 computer system at the Langley Research Center was used 
in the real-time mode for simulating the vehicle dynamics and for the on-line and off-line 
data processing. 

The feasibility of the tunnel display for operation in actual flight was demonstrated 
at NASA Wallops Flight Center, Wallops Island, Va. A modified CH-47 helicopter served 
as the flight vehicle. The tunnel image provided sufficient control information for fol- 
lowing the trajectory smoothly, and the deviations from the desired trajectory remained 
within reasonable limits. 



INTRODUCTION 


In a conventional helicopter approach under instrument flight rules (IFR), control 
information (such as position, rate, and attitude information) is presented to the pilot 
by a set of electromechanical flight instruments. On most instruments, the control vari- 
able is explicitly displayed as the deviation of a pointer from an indicated reference line. 
However, this information is rather symbolic and shows little resemblance with the 
visual world the pilot is used to. Furthermore, several instruments are used, and the 
pilot must derive the control information by scanning them. Thus, the conventional IFR 
approach requires a high level of pilot proficiency and imposes a high workload. In order 
to minimize the probability of human errors, the approach path generally is kept simple, 
with a limited number of straight and curved sections. 

In the advanced display concept under investigation, control information is pre- 
sented to the pilot in a format similar to the "through-the- windshield" visual field. A 
simplified computer -generated perspective image of the visual field is presented to the 
pilot as if it were a "tunnel in the sky," in which the three-dimensional (3-D) approach 
path is to be followed. In contrast with conventional instrument displays, the control 
information is integrated into one format, natural to the pilot. This characteristic makes 
the display very suitable for following, with great accuracy, approaches which are com- 
plicated, strongly curved, and steep. This type of approach might become a necessity in 
heavily congested fixed-wing and helicopter traffic areas. 

Use of trade names or names of manufacturers in this report does not constitute an 
official endorsement of such products or manufacturers, either expressed or implied, by 
the National Aeronautics and Space Administration. 

SYMBOLS AND ABBREVIATIONS 
A^Ay lateral or vertical not-augmented system matrix 

A^Ay lateral or vertical augmented system matrix, SAS engaged 

inertial accelerations of vehicle along y^-axis and z^-axis, respectively, 
m/sec^ 

centripetal trajectory acceleration, m/sec^ 

B^,B y lateral or vertical input matrix 
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autocovar iance 
predictor distance, m 

lateral or vertical deviation from trajectory, m 
preview term, m, see equation (A34) 

lateral intercept angle of velocity vector and trajectory direction, rad 
vertical intercept angle of vehicle axis and trajectory direction, rad 
lateral intercept angle of vehicle axis and trajectory direction, rad 
lateral or vertical feedback matrix 
lateral trajectory curvature, 

roll rate on cyclic SAS gain and feedback, respectively 

roll rate on differential, cyclic SAS gain and feedback, respectively 

pitch rate on SAS gain and feedback, respectively 

yaw rate on SAS gain and feedback, respectively 

perturbed forward-velocity SAS gain and feedback, respectively 

sideslip SAS gain and feedback, respectively 

pitch SAS gain and feedback, respectively 

acceleration due to gravity, m/sec^ 


H> 


i £ (s),H fi v (s) 


transfer function relating deviation to control input 


vehicle moments of inertia about x^-axis, y^-axis, and z^-axis, 
respectively, kg-m^ 


*xx’ 



Ij^. vehicle product of inertia in x^z^ -plane, kg-m^ 

Kj lateral gain, stability augmentation system 

d £ 

K- vertical gain, stability augmentation system 

a v 

^M’^bv gain of biases on estimates of v and w, respectively, see equations (7) 
and (8 ) 

kj 3r undesired rudder input gain, see equation (9) 

kg,k v lateral or vertical automatic control system gain, see equations (3) and (4) 

kgp simulated roll-rate feedback failure gain, see equation (10) 

kp v gain for reduction of vertical acceleration in computing vertical prediction, 

see equation (B5) 

L,M,N total rolling, pitching, and yawing moments about center of gravity of vehicle, 
N • m 

m vehicle mass, kg 

p,q,r vehicle angular velocity about x^-axis, y^-axis, and z^-axis, respectively, 

rad/sec 

R^,R y radius of curvature of lateral or vertical vehicle path, m 

Rj. radius of curvature of lateral trajectory, m 

s Laplace operator 

sc( ) deviation score for entry maneuvers, m 

T c time elapsed between moment of manual takeover and moment of reengage- 

ment of automatic control system (correction time), sec 

Tj time of occurrence of automatic control failure, sec 

T prediction time, sec 

Jr 
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T r time elapsed between occurrence of automatic control failure and moment of 

manual takeover (reaction time), sec 

T sr T sv lateral or vertical settling time, sec 

starting point and end point, respectively, of run-time range (used for 
computation of deviation scores) 

t time 

u,v,w velocity of vehicle along x^-axis, y^-axis, and z^-axis, respectively, m/sec 

U N ,W N nominal trim values of u and w, m/sec 

V ground speed, m/sec 

V vehicle velocity vector 

v ,w lateral or vertical gust component, m/sec 

O O 

W r W y lateral or vertical disturbance matrix 

X,Y, Z total aerodynamic, thrust, and gravitational forces along x^-axis, y^-axis, 
and z^-axis, respectively, N 

x b ,y b ,Zb orthogonal body-axis coordinate system, with xj^-axis pointing toward front 
of vehicle, y b -axis pointing to right, and z^-axis pointing downward 

Xj,y^,z^ orthogonal inertial coordinate system, with x^-axis pointing to north, y b -axis 
to east, and z b -axis downward 

lateral or vertical state vector 

^b’^c’^s’^r differential-collective, collective, cyclic, and differential-cyclic control 
inputs, cm 

6 b,a’ 6 c,a’ 6 s,a’^r a differential-collective, collective, cyclic, and differential-cyclic 

actuator inputs, cm 
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5^,5 y lateral or vertical control input vector 

6 0 „A, - lateral or vertical actuator input vector 

S, a V, U 

e^, e y lateral or vertical predicted deviation from trajectory, m 

£ damping factor 

0,0, 0 pitch, roll, and yaw attitude angles, rad 

6 N nominal trim value of pitch attitude, rad 

0^ downslope of trajectory, rad 

lateral or vertical displacement of predictor on screen, in screen coordinates 
cr standard deviation (root-mean-square deviation from mean) 

natural frequency, rad/sec 
Abbreviations: 

EADI electronic attitude director indicator 

FDRS flight display research system 

HSI horizontal situation indicator 

IFR instrument flight rules 

SAS stability augmentation system 

VAC vehicle axis cross 

VALT VTOL approach and landing technology 

An arrow over a symbol denotes a vector quantity. A dot over a symbol denotes the 
first derivative with respect to time. A prime denotes values generated by the automatic 
control system. 
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DISPLAY CONFIGURATION 


The conceptual tunnel display is shown in figure 1. The trajectory to be followed is 
presented as a perspective image of a tunnel with a square cross section. The perspective 
appearance of the tunnel changes with the vehicle position. Positional cues are derived 
from this perspective appearance, while rate cues are derived from the rate of change in 
the appearance. 

The tunnel cross section is constant and square with a scale width of 45.7 m (150 ft), 
which corresponds to the width of a typical runway. The Adage Graphics Terminal is 
programmed to draw (see fig. 2(a)) tunnel squares (l) and cornerlines (2) . An experi- 
mental configuration was chosen in which the tunnel squares are drawn a scale distance 
of 91.4 m (300 ft) apart. In order to augment the impression of forward velocity, the 
cornerlines are dashed. Only the five squares nearest to the vehicle are drawn. Thus, 
the number of lines drawn is substantially reduced, which prevents cluttering of the dis- 
play. The tunnel is shown in the viewing range from 0 to 762 m (2500 ft). The maximum 
horizontal and vertical field of view is from +45° to -45°. The tunnel image generated 
is square and is shown on the 33-cm (13-in.) display monitor, with aspect ratio 4:3 
(horizontal: vertical). The tunnel image magnification is 0.165 for an average viewing 
distance of 76.2 cm (30 in.) from the monitor. Pitch scales (j3) and pitch pointers (IT) 
are located at the left and right sides of the image, and a roll scale ( 5 ) is located at the 
top. Pitch motion is also visualized by an equivalent rotation of the optical axis. For 
small angles, this is equivalent to a vertical displacement of 0.22 cm per degree pitch 
angle of the tunnel image. The pitch pointer and the tunnel image are displaced, while 
the pitch scales remain vertically fixed. 

Digital readouts of magnetic heading (MAG) (To) and airspeed (TAS) (Tl) are dis- 
played in boxes on the vertical centerline at the top and bottom of the image, respectively. 
(See fig. 2(a).) These boxes also provide a reference for estimating the center of the 
image. Note that a large area in the center of the display is kept clear of symbology to 
prevent cluttering the tunnel image. 

In addition to che basic tunnel display, several other displays were used in the evalua- 
tion, either as additions to the tunnel display to improve pilot performance or as primary 
displays for purposes of comparison. The following display configurations were 
investigated. 


Roll-Stabilized Tunnel (Configuration A) 

In visualizing the roll motion, two versions were considered. In the first version 
(fig. 2(a)), the tunnel image, pitch scales, horizon, and roll scales are roll- stabilized, 
and- the bank angle <fi is visualized by rotating the wing bars (7) and the roll pointer ( 9 ) , 
in fixed-wing aircraft. This image does not conform to the through- the -windshield 
visual field and can only be used in head-down type displays. 
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Rolling Tunnel (Configuration B) 

In the second version (see fig. 2(b)), the tunnel image is rotated with the pitch and 
roll scales, and the aircraft reference symbology remains stationary. Note that, in this 
configuration, the pitch motion results in an image displacement perpendicular to the 
rotated horizon represented by the pitch pointers. 

Predictor /Director Display (Configuration C) 

In order to obtain the rate information required, it is necessary to estimate the pre- 
dicted vehicle path. As discussed in references 1 and 2, peripheral vision is utilized in 
estimating the vehicle path from the apparent motion in the visual field. Since peripheral 
vision is partially missing due to the limited size of the monitor, this estimation process 
is expected to be limited. A possible solution for this problem is to display the pre- 
dicted vehicle position explicitly on the tunnel image. A display was designed in which a 
flight-path predictor and a flight-path director are displayed. (See fig. 2(c).) The pre- 
dictor is displayed as a cross located on the predicted flight path at a distance D ahead 
of the vehicle. The director is displayed as a square with its center located on the 
desired flight path a distance D ahead of the vehicle. The outer dimensions of the 
cross and the square are equal (i.e., width: 45.7 m (150 ft)). The distance D can be 

varied, and the apparent size of the cross and square on the image varies accordingly. 

The predictor cross remains parallel to the vehicle wing line at all times, and thus, will 
remain upright in the roll version and rotate in the roll-stabilized version. With this 
display, the control task reduces to a simple tracking task. 

In a well-defined trajectory-following task, in the presence of random disturbances, 
and for a well chosen prediction time Tp, this predictor /director (i.e., predictor cross 
plus director square) display is expected to perform equally well as a more complex 
tunnel plus predictor /director (i.e., tunnel plus predictor cross plus director square) 
display. However, the predictor /director display yields the following shortcomings: 

(1) Without the tunnel image, the pilot has no knowledge of the position or attitude of 
the vehicle with respect to the tunnel. Thus, the pilot has no way of judging how well the 
predictor /director performs in following the trajectory. 

(2) With the predictor /director display, the pilot is forced to use the control error 
between the cross and the square. With the tunnel plus predictor /director display, the 
pilot has the option of partially or fully rejecting the information provided by the 
predictor /director and relying more on the tunnel image. 

(3) It is clear from shortcoming (2) that the predictor /director display will perform 
adequately only in the well-defined stationary control task for which Tp is adjusted cor- 
rectly. However, the control task generally is not stationary and is subject to large vari- 
ations in velocity, disturbance spectrum, stability derivatives, etc. 
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Tunnel Plus Predictor /Director (Configuration D) 

This configuration (see fig. 2(d)) is a combination of the roll -stabilized tunnel 
(configuration A) and a superimposed predictor /director (configuration B). The predictor 
cross ^2) shows the predicted vehicle position at a distance D ahead of the vehicle. 

The bright director square (15) represents the desired future vehicle position at dis- 
tance D. In contrast with the squares of the tunnel, this square moves along with the 
vehicle at a fixed distance D ahead. 

Vehicle Axis Cross (VAC) /Director (Configuration E) 

In this display (see fig. 2(e)), the cross represents the vehicle axis at a distance D 
ahead. The cross is locked in the center of the image and rolls with the vehicle. The 
director square is identical to the director square of configuration C. 

Tunnel Plus Vehicle Axis Cross/Director (Configuration F) 

This configuration (see fig. 2(f)) is a superposition of the tunnel (configuration A) 
and the VAC/director (configuration E). 

Combination Vertical and Horizontal Situation (CVHS) Displays 

(Configuration G) 

In this configuration (see figs. 2(g) and 2(h)), the information of conventional elec- 
tromechanical instruments is presented on two monitors. The configuration is not 
intended as a new development, but serves rather as a baseline for comparing conven- 
tional with more advanced display concepts. The upper monitor (fig. 2(g)) shows the 
horizon bar (3), pitch scales (D, roll scale ©, fixed vehicle axis reference cross 0, 
wing line ©, vertical speed indicator (§), actual altitude (2), commanded altitude d), 
and digital read outs of magnetic heading d) and airspeed (10) . In addition to this, the 
actual lateral and vertical deviation is visualized by the square (ll) , which is of the same 
size as the vehicle reference cross (§) and which corresponds to the apparent size of the 
tunnel cross section 137 m (450 ft) ahead of the vehicle. A displacement of the square, 
equal to the width of the square, corresponds to a deviation of 45.7 m (150 ft). 

On the lower monitor (fig. 2(h)), the horizontal situation is shown with vehicle ref- 
erence symbol (1^), which remains stationary, and a moving map (1J), which shows a plan 
view of the desired trajectory. 
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VEHICLE MODEL 

A linear time -invariant vehicle model, derived from reference 3, was used in the 
calculation of vehicle motion. The longitudinal and lateral dynamics and kinematics 
were assumed to be fully decoupled and linearized. The stability derivatives used in this 
model are valid for relatively small perturbations about a trim condition in which the 
vehicle is in a level flight with a commanded null velocity of V = 33.5 m/sec (110 ft/sec). 
A detailed description of the vehicle model and the state equations for the nonaugmented 
lateral and longitudinal dynamics are given in appendix A. The lateral and longitudinal 
stability derivatives are given in table 1. 

EXPERIMENTAL EVALUATION 

Objectives of the Experimental Program 

The objectives of the experimental program are: 

(1) To determine whether the three-dimensional tunnel image provides sufficient 

positional, rate, and attitude information for entering and following a steep 

and curved trajectory in a moderate to heavy gust environment 

(2) To investigate the effect of predictive information 

(3) To compare the performance of the tunnel display with the performance of 

displays based on conventional approach instrumentation 

(4) To investigate the effect of roll cues 

(5) To investigate the abilities of the tunnel display in monitoring automatic 

approaches 


Description of the Experimental System 

The experimental program was carried out on the Langley Real-Time Simulation 
System. A functional diagram of the experimental system is shown in figure 3. A 
Control Data CYBER 175 computer system at the Langley Research Center was used on a 
time-sharing basis for the real-time digital simulation of the vehicle dynamics. The 
computed vehicle motions were imparted to an Adage Graphics Terminal, the AGT/130, 
which was programmed to generate images of the various display configurations. These 
images were displayed on two display monitors, mounted one above the other, in a fixed- 
base simulator cabin equipped with helicopter controls (fig. 4). Vehicle motions, thus 
presented on the cabin display monitors, were utilized by the pilot to generate the control 
commands, which were imparted to the CYBER computer system. 
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Prior to each experimental run, the experimental conditions, such as display 
parameters, run control parameters, initial conditions, and disturbances, were set at 
the CYBER 175 console. On-line data processing included the computation of actual and 
predicted deviations from the trajectory as well as the computation of running averages 
of these deviations, of the vehicle motions, and of the control commands. 

Note that all display computations have been performed by the Adage Graphics 
Terminal on a stand-alone basis, without interfacing with the CYBER 175 computer. The 
only variables transferred from the AGT/130 to the CYBER 175 were the vehicle motions 
and the control parameters which defined the display to be generated. 

The reasons for operating the AGT/130 on a stand-alone basis are twofold: 

(1) It enabled direct implementation of the display in the Flight Display Research 
System (FDRS) at the NASA Wallops Flight Center (ref. 4), where a real flight vehicle was 
used and measurements of actual vehicle motions were fed to the AGT/130. 

(2) Stand-alone graphics techniques, thus developed, prepared the way for the 
development of techniques for onboard generation of the display with special purpose 
microprocessors. 

The real-time simulation program was written in Fortran and used 102g K memory 
storage. A fourth-order Runge-Kutta integration scheme was used for solving the equa- 
tions of motion. Integration, sampling, and averaging was at a rate of 32 Hz. Data could 
be recorded at 32 Hz as well. Overall computation time of one iteration amounted to 
40 percent of the time interval. Data were transferred from the CYBER 175 to the 
AGT/130 at a rate of 32 Hz in 15-bit fixed-point format. The images generated by the 
AGT/130 were refreshed at 30 Hz, but the tunnel image was updated at 12 Hz. 

The simulator cabin was equipped with a two-axis electrohydraulic control stick, 
spring loaded rudder pedals, and an unloaded, balanced, collective control lever. The 
range of the two-axis control stick was from +12.7 to -12.7 cm (+5 to -5 in.), and the 
spring torque was 0.0848 N • m (0.75 lbf • in.) for both axes. The range of the collective 
level was from +12.7 to -12.7 cm (+5 to -5 in.). As discussed previously, the zero-stick- 
force velocity of the vehicle was kept constant at 33.5 m/sec (110 ft/sec), and the forward 
and aft control stick motion generated a command to increase or decrease the velocity 
about 0.75 m/sec per cm stick deflection (6 ft/sec per in. stick deflection). The lateral 
control stick motion generated a roll-rate command, resulting in an initial roll rate of 
about 0.058 rad /sec per cm stick deflection (8.5 deg/sec per in. stick deflection). Turn 
coordination was carried out automatically, and the rudder pedals were not used. The 
vertical speed was controlled by the collective lever, geared to a vertical speed at 
2.04 m/sec per cm lever deflection (17 ft/sec per in. lever deflection). Null position of 
the collective lever corresponded to zero vertical velocity. 
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The electromechanical instruments of figure 4 were not connected, and all control 
information was derived from the two display monitors. On the upper monitor (measuring 
33 cm (13 in.) diagonally), a raster-scan image of the tunnel, or of an electronic attitude 
director indicator (EADI), was shown. On the lower monitor (measuring 22.9 cm (9 in.) 
diagonally), a horizontal situation indicator was shown in part of the experiments. The 
average distance from the subjects eye to the monitors was about 76.2 cm (30 in.). 


Description of the Experiments 

All experiments were concerned with the approach phase of the flight in the range 
from 4570 to 152 m (15 000 to 500 ft) from the landing point. A plan view of the desired 
trajectory is shown in figure 5, and the desired vertical descent profile along the tra- 
jectory is shown in figure 6. In all experiments, the forward velocity was set at 
33.5 m/sec (110 ft/sec) with vertical speed, bank angle, and pitch angle initially set 
to zero. 

The experiments were designed to cover a broad range of frequently occurring con- 
trol situations. Following are discussions of the three types of experiments that were 
chosen. 

Trajectory following in the presence of random, lateral, and vertica l gust 
disturbances .- The pilot was instructed to minimize the lateral and vertical deviations 
from the trajectory with minimum control effort. Each run started from initial location 7 
(see fig. 5) with a lateral deviation of 30.5 m (100 ft) to the right of the trajectory, a ver- 
tical deviation of zero, and an intercept angle of zero. Thus, control action was required 
from the pilot immediately after starting the run, to bring the vehicle back on the trajec- 
tory. Gust disturbances were generated by passing band-limited white noise through a 
first order filter with a break frequency of 0.2 rad/sec. The root-mean-square value of 
the lateral gusts was 5.33 m/sec and of the vertical gusts 4.57 m/sec. 

Each run lasted 120 sec, during which time the means and autocovariances of devia- 
tions, vehicle motions, and control commands were computed according to the equation 


cov(x) 


i f* 
■t Jo 


x z dt 


( 1 ) 


where cov(x) is the autocovariance of variable x, and x can be p, d^., d c , d , 


'It 9 V V 


or 5 


Entering the trajectory from a r andomly chosen, unknown location outside the 
tunnel . - This experiment attempted to simulate the situation in which the pilot, after not 
viewing the display for some time, was at once confronted with the situation of being 
located outside the tunnel. The pilot was instructed to bring the vehicle back on the 


12 



desired trajectory as fast as possible, as smoothly as possible, and with minimum con- 
trol effort. Gust disturbances were not present in this experiment. In order to prevent 
the pilot from knowing his initial position before the start of a simulation run, the tunnel 
image, the tunnel cross-sectional square, the map, and the altimeter scale were initially 
blanked and were only made visible immediately after starting the simulation run. 

Each run started randomly from one of the six initial locations shown in figure 5. 
For all locations, the initial lateral deviation was 305 m (1000 ft) to the left or to the 
right of the trajectory, and the initial vertical deviation was 45.7 m (150 ft) above or 
below the trajectory. The initial intercept angle was set between 0° and 72°, 

Each entry run lasted 50 sec, during which the following performance scores were 
computed: the lateral settling time T g ^ - defined as the time from the start of the run 

to the moment the lateral deviation remains within a ±15.2-m (50-ft) settling tolerance 
about the desired trajectory - and the vertical settling time T gy - for which the settling 
tolerance is ±10.7 m (35 ft); the autocovariances of the roll rate p, of the intercept 
angle d^, of the predicted lateral and vertical deviations and e y and their rates 
and e y , and of the control commands 6 g and 6 C ; the lateral deviation score - 
defined as the averaged absolute value of the lateral deviation computed between t = 
and t = T 2 , where = 10 sec and T 2 = 50 sec according to 

° c w*(^C>i di <2> 

The vertical deviation score was computed in the same way, with T^ = 6 sec and 
T 2 = 50 sec. The lateral and vertical deviation scores were chosen to be averaged abso- 
lute values rather than averaged squared values, in order to prevent these scores from 
being dominated by the large initial deviation. For the same reason, the averaging pro- 
cess only starts at t = T^ sec. Starting point T^ was chosen to be 25 percent less 
than the best previously recorded settling time. 

Monitoring experiment . - The purpose of this experiment was to investigate the 
abilities of a pilot, using the various display configurations, to detect and correct control 
system failures in automatic approaches. Each run started as an automatic flight. At a 
random chosen time t = Tj (between 5 and 20 sec), a failure occurred in the automatic 
control system. After detecting the failure, the pilot responded by pressing one of two 
buttons on the control stick. This disengaged the automatic control system and engaged 
the manual control mode. The control mode, automatic or manual, was indicated on the 
display by AUTO or MAN, respectively. The pilot was instructed to bring the vehicle 
back on the desired trajectory as fast as possible and with minimum control activity. 

After returning the vehicle to the desired trajectory, the pilot had to reengage the 
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automatic control system by pressing the second button on the control stick. The control 
failure was assumed to be eliminated by the use of a backup automatic control system. 
The condition for the successful reengagement of the automatic control system was that 
the lateral and vertical deviations and the intercept angle must be within preset toler- 
ances. These tolerances were ±19.8 m (65 ft) for the lateral deviation, ±15.2 m (50 ft) 
for the vertical deviation, and ±0.15 rad for the intercept angle. 

Each run lasted 50 sec, during which the following performance scores were mea- 
sured: the reaction time T r - defined as the time elapsed between the occurrence of 
the failure and the moment of manual takeover; the correction time T c - defined as the 
time elapsed between the moment of manual takeover and the moment of reengagement 
of the automatic control system; the autocovariances of d £ , d^, e^, e y , d^, p, 6 g , 

and 6 p measured over the total 50 -sec interval of combined automatic and manual 
control. 

Based on the predicted lateral and vertical deviations at 229 m (750 ft) ahead of 
the vehicle, the lateral and vertical automatic control laws were 


V = -Vv < 4 > 

where automatic control inputs 6 g ' and 6 ' replaced the manual control inputs 6 g 
and 6 C to the system with SAS engaged. The gains of the automatic control system were 
set to 0.018 for k^ and 0.045 for ky, and the vertical predictor gain kp V was 0.2. 

During each run, a control failure would occur, chosen at random from the four 
failures described hereafter: 

(1) Lock of control inputs, where 


V<«=V( T f) 

(5) 

V<‘> =V( T f) 

(6) 


for t£ T f . 

(2) Error estimates for v and w, which were used in computing the predicted 
vehicle position. In this study v and w were assumed to be perfect estimates. How- 
ever, a control failure would occur when incorrect values of v and w were used in 
computing the predicted vehicle position. Additive errors, linearly increasing with time, 
were introduced. The error estimate values v and w are then given by 
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( 1 ) 


5 - T + *bt (‘ - T f) 

" = W + k bv (t - T f ) (8) 

where k^ = k^ = 0 for t< Tj with k^ =±5 and k bv = ±2 for tS Tf. 

(3) Undesired linearly increasing rudder input, according to 

V - *br(‘ - T f) < 9 > 

where k br = 0 for t < T f and = ±0. 1 for t ^ T f . 

(4) Roll-rate feedback failure, simulated by 

V = “Vi + k £pP ( 1Q ) 

where =0 for t < T f and k £ = 7 for t = T f . 


Results 

Four subjects participated in the evaluation program- The background of the sub- 
jects is listed in table 2. Subject training consisted of two to four simulator sessions of 
1 to 3 hours duration each, in which the subjects familiarized themselves with the various 
display configurations. In addition to this, 3 to 12 training runs were made, prior to the 
production runs, for each experimental condition. The four subjects completed a total 
of about 1640 production runs, totaling 28.6 hours of net simulation time. Each subject 
performed at least four repetitions for each experimental condition. The results for the 
four subjects were treated separately and are summarized in tables 3 to 12. The result 
for each experimental condition represents the average and standard deviations of a set 
of four or more repetitions. 

Significance tests on the difference between averages were performed with a 
Student-t test at a 5 -percent level of significance. Although significant differences 
between the results of the four subjects were observed, the general trends were rather 
similar. Therefore, only the results of subject 1 were plotted. 

A 40-min motion picture was made, in which the most characteristic results of the 
three experiments are recorded. Since the motion picture was made in real time 
directly from the display monitor and during actual simulation runs, it accurately repro- 
duces the simulator experiments. 
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Results for th e traj ecto ry following experiment. - The results of the trajectory fol- 
lowing experiment for the four subjects are summarized in tables 3 to 7. The results for 
subject 1 are shown in figures 7 to 22. From these results, the following can be 
concluded. 

Comparison of the conventional EADI/map display with the basic tunnel display: 

For all four subjects, the tunnel display yielded a markedly better accuracy in following 

the trajectory than the EADI/map display, as can be seen from the significantly lower 

values for the autocovariances of d fJ and d^ for the tunnel display. (See tables 3 to 5 

and figs. 7 and 8.) The autocovariance of d . was also significantly lower for the tunnel 

V' 

display, indicating a less oscillatory lateral response. No significant difference in the 
roll activity, as indicated by cov(p), was found between the two displays. (See fig. 9.) 

These results are confirmed by the time histories of the lateral and vertical 
responses, shown in figures 10 and 11. Note the stepwise character of the collective con- 
trol input for the EADI/map display in figure 11. This might result from the fact that the 
pilot intermittently controls the lateral and vertical situations. 

Comparison of the basic tunnel with the tunnel/vehicle axis cross (VAC): The 

effect of adding a VAC and a director to the tunnel image was investigated for subjects 1 
and 2 with the cross and director at D = 305 m (1000 ft). 

No significant difference in the results was found between the basic tunnel and the 
tunnel plus VAC/director. (See tables 3 and 4 and figs. 7 to 9.) Thus, a VAC /director 
does not contribute to a better performance. 

Comparison of the tunnel plus VAC/director with the VAC/director: The 

VAC/director at D = 305 m (1000 ft) was investigated for subjects 1 and 2. The auto- 
covariances of the lateral and vertical deviations were significantly larger for the 
VAC/director than for the tunnel plus VAC/director. (See tables 3 and 4 and figs. 7 
and 8.) This means that the error between the VAC and the director does not provide 
the correct control information for this task. Thus, with the tunnel plus VAC/director, 
all control cues are derived from the tunnel image and the VAC is ignored. This assump- 
tion is confirmed by the significantly larger value for the autocovariance of for the 
tunnel plus VAC/director. 

Comparison of the basic tunnel with the tunnel plus predictor /director: For all 

four subjects, performance was found to improve significantly with the addition of 
predictor/director symbols to the basic tunnel display, as can be concluded from the 
markedly reduced roll activity and smaller lateral and vertical deviations. (See 
tables 3 to 5 and figs. 7 to 9.) 

In figures 10 and 11, the time histories for the tunnel plus predictor/director are 
compared with those for the basic tunnel. Note the lower roll and pitch activity and the 
smaller vertical deviations for the tunnel plus predictor /director. 
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The effect of the predictor distance D: The effect of the predictor distance D 

was investigated for subjects 1 and 2 for the predictor with kp y = 1 and D = 91.4 m 
(300 ft), 137 m (450 ft), 229 m (750 ft), and 305 m (1000 ft). Lateral and vertical devia- 
tions were found to increase strongly with D. (See figs. 12 and 13 and tables 6 and 7.) 

On the other hand, the roll activity decreased with increasing D. (See fig. 14 and 
tables 6 and 7.) These results indicate that at D = 91.4 m (300 ft) the system is more 
oscillatory, indicating lower damping, and at D = 305 m (1000 ft) the system is too 
inaccurate. For this task and airspeed of 33.5 m/sec (110 ft/sec), the optimum predictor 
distance tends to be between D = 137 m (450 ft) and D = 229 m (750 ft). 

Note that the autocovariance for the predicted vertical error rate e y strongly 
increases with D. (See tables 6 and 7.) This results from the fact that the vertical dis- 
turbances are less filtered by the vehicle dynamics than the horizontal ones, which 
results in rapid vertical predictor motions. 

The time histories for the tunnel plus predictor /director at D = 91.4 m (300 ft) 
and D = 229 m (750 ft) are shown in figures 15 and 16 and confirm the results mentioned 
previously. 

Comparison of the tunnel plus predictor /director with the predictor /director: The 

predictor /director was investigated for subjects 1 and 2 and for D = 91.4 m (300 ft), 

229 m (750 ft), and 305 m (1000 ft). At D = 229 m (750 ft) and 305 m (1000 ft), the auto- 
covariances of the lateral and vertical deviations for the tunnel plus predictor /director 
were found to be generally smaller than for the predictor/director. (See tables 6 and 7 
and figs. 12 and 13.) At D = 305 m (1000 ft), the differences were found to be significant. 
The autocovariance of the predicted lateral deviation was significantly larger with 
the tunnel plus predictor /director than with the predictor /director. (See tables 6 and 7.) 
This indicates that, for a wrongly adjusted predictor distance, the pilot relies on the 
tunnel image rather than on the error between predictor and director. Note also the sig- 
nificantly lower roll activity for the predictor/director at D = 305 m (1000 ft). 

The time histories for the tunnel plus predictor /director and for the predictor/ 
director, both at D = 305 m (1000 ft), are shown in figures 17 and 18 and confirm these 
findings. 

The simplified predictor and the effect of a reduced predictor gain: The simplified 

predictor, based on the bank angle <£, was investigated for subjects 1 and 2 with 
D = 229 m (750 ft). The autocovariance of the lateral deviation was significantly larger 
with the simplified predictor than with the regular one, but it was still significantly 
smaller than with the tunnel without predictor/director. (See tables 3 and 4.) No signif- 
icant difference in roll activity was found between the simplified predictor and the regular 
one. These results indicate that the simplified predictor yields a somewhat less accurate 
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but adequately damped system. Thus, the roll angle may very well be used for the lateral 
prediction, in those situations where the lateral acceleration is not available. 

In the vertical control, the simplified predictor with kp V = 0 yielded vertical 
deviations similar to those of the predictor with kp V =0.2 and with kp v = 1. However, 
because of the rapid vertical predictor motions, the rate of the vertical error between the 
predictor and the director e v was found to be very large for kp V = 1. The test subjects 
objected to these rapid predictor motions. On the other hand, the vertical accelerations 
were completely disregarded for kp V = 0, which resulted in overcontrolling because of 
the lack of quickening in the display. A compromise between sufficient quickening and 
acceptable vertical predictor motions was found for kp V = 0.2. 

The effect of roll cues: The effect of roll cues was investigated for subjects 1 

and 2 with the tunnel plus predictor /director at D = 229 m (750 ft) and kp y = 0.2. 

No significant difference was found in the results between the roll version and the roll- 
stabilized version. (See tables 3 and 4.) This indicates that, for this control task, the 
roll version and the roll-stabilized version perform equally well in a fixed-base 
simulator. 

Results for the trajectory en try experiments. - The entry experiments were con- 
ducted in a series of six runs. For each run, the initial location was chosen at random 
from one of the six positions given in figure 5. No position was used more than once. 
Significant differences between the scores were observed for the various initial locations. 
In order to rate the general entry performance, the results of the six runs in each series 
were averaged in order to obtain the series scores. For each display configuration, two 
to six series of runs were performed. The experimental results for the four subjects 
are summarized in tables 8 to 11. Each result represents the average and standard 
deviation of a set of four to six series scores. From these results, the following can be 
concluded. 

Comparison of the EADI/map display with the basic tunnel display (roll version): 

For all four subjects, the entry performance for the basic tunnel was found to be markedly 
better than for the conventional EADI/map display. With the basic tunnel, significantly 
smaller lateral and vertical deviation scores, as well as smaller lateral and vertical 
settling times, were found (see tables 8 to 11 and figs. 19 to 22), indicating a faster and 
more accurate entry. For the tunnel display, the roll activity was found to be signifi- 
cantly lower as well, indicating better system damping. (See fig. 23.) 

Shown in figures 28 and 29 are the time histories for the entry from location 3, for 
the EADI/map, and for the basic tunnel. Both displays yielded considerable roll activity, 
as shown in figure 24. However, the tunnel display yields a more gradual and accurate 
entry, as can be seen from the time history of the lateral deviation. Similar results in 
the vertical control are shown in figure 25. Note the stepwise character of the collective 
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control input, for the EADI/map. This phenomenon was discussed earlier for the tra- 
jectory following experiment. 

Comparison of the basic tunnel with the tunnel plus predictor /director: The effect 

of the predictor is noticed mainly in the markedly lower roll activity, as seen in figure 23. 
For the tunnel plus predictor /director, lateral and vertical deviation scores, as well as 
lateral and vertical settling times, were found to be significantly lower also. (See figs. 19 
to 22.) Thus, the tunnel plus predictor /director yields a more accurate and better damped 
entry. 

The effect of the predictor distance D: Lateral and vertical deviation scores, as 

well as lateral and vertical settling times, were found to increase with increasing D. 

The increase of these scores was not significant in the range between D = 91.4 m (300 ft) 
and D = 229 m (750 ft). However, a significant increase in these scores was found 
between D = 229 m (750 ft) and D = 305 m (1000 ft). This trend is shown in figures 26 
to 29. On the other hand, the roll activity was found to decrease with increasing D. (See 
fig. 30.) This decrease in roll activity was found to be significant between D = 91.4 m 
(300 ft) and D = 137 m (450 ft). These results indicate that, at D = 91.4 m (300 ft), the 
system lacks damping and, at D = 305 m (1000 ft), the system is too inaccurate. It is 
clear that, for the given forward velocity of V = 33.6 m/sec (110 ft/sec), a predictor 
distance in the range from D = 137 m (450 ft) to D = 229 m (750 ft) yields the best 
compromise. 

In figures 31 and 32, the time histories for the entry from location 2 are shown for 
the tunnel plus predictor at D = 91.4 m (300 ft), 137 m (450 ft), and 229 m (750 ft). The 
time histories of the lateral deviation and intercept angle indicate a smoother entry for 
larger D. (See fig. 31.) The time histories of the roll rate for D = 91.4 m (300 ft) 
indicate a considerable increase in roll activity at the point of intercept with the tunnel. 
The lack of damping at D = 91.4 m (300 ft) was less pronounced in the vertical control. 

It was noticed in a slightly increased pitch rate. (See fig. 32.) 

Comparison of the tunnel plus predictor /director with the predictor /director: The 

predictor /director display was investigated for D = 229 m (750 ft) and 305 m (1000 ft). 
Predictor distances smaller than D = 229 m (750 ft) were not attempted, since the 
tunnel square was not visible during the first part of the entry. 

At D = 229 m (750 ft) and 305 m (1000 ft), lateral and vertical deviation scores, 
as well as lateral and vertical settling times, were found to be significantly larger for the 
predictor /director than for the tunnel plus predictor /director. (See figs. 26 to 29.) This 
proves that the tunnel image contributes significantly to the entry accuracy, particularly 
for a badly adjusted predictor (e.g., D = 305 m (1000 ft)). The roll activity for the tunnel 
plus predictor /director was somewhat larger than for the predictor/director. When the 
tunnel image was shown in addition to the predictor /director, the pilot relied on the tunnel 
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image rather than on the error between predictor and director square. This is seen in 
the fact that the autocovariance of the predicted lateral error was significantly larger 
for tunnel plus predictor /director than for predictor /director. 

In figures 33 and 34, the time histories for the entry from location 4 are shown for 
the tunnel plus predictor /director and for the predictor /director. It is shown in figure 33 
that, for the tunnel plus predictor/director, the pilot allows a rather large predicted 
lateral error, while the time history of the actual lateral deviation shows a smooth and 
accurate response. On the other hand, the predicted lateral error was kept small for the 
predictor /director; whereas, the time histories of the actual lateral deviation show a 
large overshoot. 

Similar results were obtained for the vertical control, as can be seen from the 
time histories in figure 34. 

The effect of roll cues: The roll and roll stabilized versions of the tunnel plus 

predictor /director at D = 137 m (450 ft) were compared, and the results are shown in 
tables 8 and 9. The differences between the scores were not significant, even at a 
10-percent confidence level. This indicates that, in the entry task, the roll and roll- 
stabilized versions perform equally well for fixed-base simulation. 

Results for the monitorin g ejcperiment. - The monitoring experiments were con- 
ducted in a series of four runs of 50 sec each. In each run, one control failure appeared, 
chosen at random from the set of four failures discussed earlier. The polarity of the 
failure biases were chosen at random as well. The results of four runs in a series were 
averaged, in order to obtain averaged failure scores for the series (series scores). Sets 
of at least 12 series of four runs were conducted, for each experimental condition. 

The monitoring experiment included the tunnel plus predictor/director at D = 229 m 
(750 ft), the predictor /director at D = 229 m, and the EADI/map. The results for sub- 
jects 1 and 2 are summarized in table 12. Each result represents the average and stand- 
ard deviation of a set of 12 to 19 series scores. These results show that, for the lateral 
deviation, the vertical deviation, as well as the intercept angle, the tunnel plus predictor/ 
director yielded significantly smaller overall autocovariances than the other two displays. 
The EADI/map performed worst in this respect. The same trend is noticeable for the 
reaction time T r and the correction time T c . The values of T f and T c were 
found to be significantly smaller for the tunnel plus predictor/director than for the 
EADI/map. Thus, the tunnel plus predictor /director enabled the best failure discovery 
and failure correction. 


Qualitative Flight Test Validation 

The feasibility of the tunnel display for operation in actual flight was demonstrated 
at the NASA Wallops Flight Center. A CH-47C helicopter (fig. 35) served as the flight 
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vehicle. A special version of a tunnel/map display was developed for operation at 
Wallops. The displays were generated by a ground-based Adage Graphics Terminal, 
scan converted, transmitted by TV link to the vehicle, and displayed on two monitors 
mounted in the CH-47C cockpit. On the upper display, the roll version of the tunnel plus 
predictor /director was shown; on the lower display, an HSI combined with an altimeter 
and a vertical speed scale. This configuration, which was specially designed for the 
VTOL approach and landing technology (VALT) program, served primarily as a display 
for conducting and monitoring automatic control system experiments. 

However, to validate the tunnel display, a limited number of manual approaches 
were conducted. Subject 3 was the research pilot, and he used the same approach trajec- 
tory as used in the simulator experiments. The CH-47C was flown with the SAS engaged, 
but the velocity hold system used in the simulator was not available. Although the 
approach trajectory was the same as in the simulator experiments, the control task was 
considerably more demanding, since the pilot also had to keep the forward velocity con- 
stant. Furthermore, the predictor /director was less effective than in the simulator 
experiments, since not all variables necessary to drive the predictor /director were avail- 
able. The lateral prediction was based on the bank angle <fi according to equation (B6). 
(See appendix B.) Since v was not available, it was set to zero in equation (B3). (See 
appendix B.) In the vertical prediction, kp V was set to zero. 

In spite of the task difficulty and the lack of training runs, the trajectory was fol- 
lowed smoothly, while the lateral and vertical deviations from the trajectory remained 
within ±30.5 m (100 ft) and 15.2 m (50 ft), respectively. 

The. pilot commented that, due to the lack of a velocity hold system, the workload 
was rather high. However, the tunnel image provided sufficient information for following 
the trajectory. Since the predictor /director did not provide the correct control cues, the 
pilot chose to utilize the tunnel image rather than the predictor /director. However, the 
predictor /director did provide the damping cues to some extent. The pilot also com- 
mented that, in spite of the task difficulty, the display was acceptable and that solving the 
problems mentioned would reduce the workload considerably. 

CONCLUSIONS 

The tunnel display with a superimposed predictor/director symbol was found to be 
a very successful combination. It outperformed the other two displays in all three exper- 
iments. A prediction time of 4 to 7 sec was found to optimize trajectory tracking for the 
given vehicle dynamics and flight conditions. Pilot acceptance of the tunnel plus predictor/ 
director display was favorable, and the time for familiarization was relatively short. 
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Specifics: 

1. The perspective tunnel image provides adequate positional and directional infor- 
mation and yields markedly better accuracy in trajectory following and trajectory entry 
than the conventional EADI/map display. 

2. The basic tunnel display yields poor damping, due to the lack of rate information 
(caused by the narrow visual field). 

3. Predictor/director symbols superimposed on the tunnel image provide the rate 
information which is vital for an adequately damped system. For the given vehicle 
dynamics and nominal forward velocity of 33.5 m/sec (110 ft/sec), a predictor symbol 
predicting the vehicle position 4 to 7 sec in advance yields the best compromise between 
positional accuracy and system damping. 

4. The predictor /director display performs well in trajectory following. However, 
the use of this display is limited to the narrow range about the trajectory in which the 
director is visible. 

5. The tunnel plus predictor /director display outperforms the predictor/director 
display when the predictor is badly adjusted, or in situations in which the director is out 
of the visible range. 

6. A prediction time of 4 to 7 sec, which yields the best compromise for the lateral 
control, yields vertical predictor motions which are too rapid. A vertical prediction 
based on 20 percent of the actual vertical acceleration yields the best pilot acceptance. 

7. A vehicle axis cross, superimposed on the tunnel image does not contribute to 
better tracking. The error between vehicle axis cross and director square does not pro- 
vide the correct control cues in trajectory following in a fixed-base simulator. 

8. The roll and roll-stabilized versions of the tunnel display perform equally well 
in trajectory following and entries. This leaves the option open, for applying the roll 
version to head-up displays, in which the image has to be conformal with the visual 
world. 

9. The tunnel plus predictor /director display performs very satisfactorily in the 
detection and correction of failures in automatic flight. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 1, 1980 
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VEHICLE MODEL 

Vehicle Model and Stability Augmentation System 

All equations of motion comply with the north-east-down sign convention. The 
lateral equations of motion, for small deviations from the nominal trim condition are 
given by 
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where 
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For small angles of 6 , the horizontal coordinates of the vehicle location in the inertial 
reference system, x i and y p are obtained by integration of 
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The longitudinal equations of motion, for small deviations from the nominal trim condi- 
tion, are given by 
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*v = \S, + B v 6 va + W v w g 


(A6) 


The altitude z i is obtained by integration of 
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Zi = V sin - flj (A7) 

A block diagram of the lateral SAS is shown in figure 36(a). The sideslip feedback is 
given by 


Sv = f v v 


(A8) 


The yaw-rate feedback is given by 
Sr = f r( r " Si) 


(A9) 


gj = -0.312gj + 0.312r 


(A10) 


The roll rate on differential cyclic (rudder) feedback is given by 

g pr = g 2 (AH) 

g 2 = -0.243g 2 + 0.243f pr p (A12) 


and the roll rate on cyclic feedback is given by 


g 


P 



(A13) 


With the SAS engaged, the cyclic and differential-cyclic actuator inputs _ and 6 

s, a 

are given by 


s,a 


6 - g 

s & p 


(A14) 


6 r,a =6 r + Spr-Sr + Sv ( A15 ) 

where 6 g and $ r are the lateral-control stick and rudder -pedal commands, 
respectively. 

State equation (Al) can be augmented, in order to incorporate the SAS dynamics. 
The lateral, augmented open-loop (SAS disengaged) state equation becomes 
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or written shortly as 


*£ = A 'A + B i&,a + W i V g 


The lateral, augmented closed-loop (SAS engaged) system is given by 


*i ' A i T i + B i 6 t * w i v g 


where A^ is the lateral closed-loop system matrix given by 



and is the lateral feedback matrix given by 


(A16) 


(A17) 


(A18) 


(A19) 
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The lateral SAS can be easily altered or disengaged by altering or nullifying the 
matrix . 

A block diagram of the longitudinal SAS is shown in figure 36(b). It consists of a 
pitch-rate feedback 

g q = f q Q (A22) 

a pitch-angle feedback 

g e =v (A 23 ) 

and a perturbed forward-velocity feedback 

gu = f u u (A24) 

The vertical closed-loop system is given by 

\ = K\ + Vv + W v w g (A25 ) 

where is the closed-loop system matrix given by 

\ - 0V - B v F^J (A26) 

and F y is the vertical feedback matrix, given by 
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0 f q 

0 0 0 


(A27) 


It is clear from the lateral stability derivatives that the not-augmented vehicle had 
very poor directional stability and poor yaw and roll damping. The lateral stability aug- 
mentation system (SAS) was derived from reference 5 and is shown in figure 36(a). It 
consists of a sideslip feedback, which provides the directional stability; a yaw-rate feed- 
back, which provides the necessary yaw damping; a roll rate on differential cyclic (rudder) 
feedback, which minimizes the effect of roll-induced yaw; and a roll rate (on cyclic) feed- 
back, to provide the necessary roll damping. With the SAS engaged, lateral control is 
accomplished by a lateral stick command only, without using the rudders, while the turn 
coordination is accomplished automatically by the feedback on the sideslip. 

The not-augmented longitudinal vehicle dynamics had an unstable pole at 
s = 0.519 rad/sec. The longitudinal SAS was specially designed for this study and is 
shown in figure 36(b). The system had to provide longitudinal stability and had to keep 
the forward velocity constant. This was accomplished by a pitch-rate, pitch-attitude, 
and perturbed-forward-velocity feedback on the differential collective input (pitch). Thus, 
a forward stick command resulted in a proportional increase or decrease of the forward 
velocity. Vertical control was accomplished by a collective command. The rotor speed 
was assumed to remain constant, and the stability derivatives were assumed to be invar- 
iant with the engine torque. 

d, 

The pole-zero location of the transfer function (s), relating the linearized 

lateral vehicle displacement with a cyclic control input, is shown in figure 37(a) for the 
SAS disengaged and in figure 37(b) for the SAS engaged. Shown in figures 38(a) and (b) 

d 

are the locations for the IT v (s) transfer function, relating the linearized vertical dis- 

0 c 

placement with a collective input. As shown in figure 37(b), the lateral vehicle dynamics 
with SAS engaged were dominated by the three poles at the origin. The effect of the pair 
of complex poles with a natural frequency o> n = 0.32 rad /sec and a damping factor of 
? = 0.43 is less dominant, because of the closely located pair of complex zeros, with 
co Q = 0.25 rad/sec and £ = 0.43. 

The longitudinal dynamics with SAS engaged, shown in figure 38(b), were dominated 
by one pole at the origin and one at s = -0.462 rad/sec. The effect of the well damped 
pair of complex poles, with a; n = 0.73 rad/sec and £ = 0.67, was masked by the closely 
located pair of complex zeros, with = 0.63 rad/sec and £ = 0.76. 
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Control Task Definition and Candidate Loop Closures 

The control task was to minimize the lateral and vertical deviations from a speci- 
fied reference trajectory, with minimum control effort. 

The horizontal situation for following an arbitrarily curved trajectory is shown in 
figure 39. The lateral deviation from the trajectory is denoted by d^, the intercept 
angle of velocity vector and trajectory direction by d^, and the curvature of the trajec- 
tory at the approximate vehicle location by f(t) = l/R t (t). The centripetal accelera- 
tion a^ for following a curve with radius R^ is given by 

2 

a t (t) = ir = V 2 f(t) (A28) 

R t w 

where V is the ground speed. For small deviations, d^ and d^ are given by 

<yt) = i J [a £ (t) - a t (t)j dt (A29) 

d £ (t) = V I d y (t) dt (A30) 

where a ^ is the acceleration of the vehicle in lateral direction and with respect to an 
inertial reference system. 

d c 

It has been shown in figure 37(b) that the transfer function tL*(s) of the lateral 

°s 

augmented vehicle dynamics was dominated by three poles at the origin. It is clear that 
adequate rate information is required in order to stabilize the system. It is shown in 
references 1 and 2 that rate cues can originate from the visual field by deriving the con- 
trol error a distance D ahead of the vehicle. It is also shown that the pilot was able to 
obtain a rough estimate of the predicted vehicle path from the apparent motion of con- 
spicuous objects in the visual field (texture points). In figure 39 this circular, predicted 
vehicle path is shown. Point A is the predicted vehicle position at a distance D, or 
Tp = D/V sec, ahead of the vehicle. The control error is defined as the deviation of 
point A from the reference trajectory. It follows from the geometry of figure 39 that, 
for small d^d^ and D^R^., the following linearized expression for can be derived: 

e £ = d £ + Dd y + |(a £ - a t ) + d p (A31) 
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or, with equations (2) and (3) and after Laplace transformation, 

e je (s) = + T p s + 1 T p 2 s 2 ^ d £ (s) + d p (s) (A32) 

The term d p constitutes a preview of the trajectory and accounts for changes in curva- 
ture d p (t) given by 

C7=Tp 

d p (t) = JJ [a t (t) - a t (cr)J d 2 CT (A33) 

0 

It follows from equation (A33) that, for straight trajectories or trajectories with a con- 
stant curvature, dp will be zero. 

A block diagram of this linearized system is shown in figure 40(a). The curva- 
ture f(t) is the forcing function of the system. 

A control error, based on the predicted vehicle position, yields the following 
advantages: 

(1) A pair of complex zeros with a fixed damping factor of £ = 0.707 and a natural 
frequency u> n which can be adjusted to the requirements of the system by varying the 
prediction time Tp. 

(2) For straight trajectories or trajectories with a constant curvature, a zero 
value of will yield a zero steady-state value of d^, as can be seen by setting e^, 
s, and dp to zero in equation (A32). 

An alternative loop closure can be formulated by considering a point B on the 
vehicle direction of motion Tp seconds ahead of the vehicle. (See fig. 39.) The con- 
trol error is then defined as the deviation of point B from the trajectory and, for small 
deviations, is given by 

e £ (t) = d £ (t) + Ddy (t) + | T p 2 a t + dp(t) (A34) 

or 

e £ (s) = ( ■* + Tp s ) df + 2 T p 2 a t (s) (A35) 

A block diagram of this system is shown in figure 40(b). Unlike the control error of 
equation (A32), equation (A35) provides only a single zero. It also follows from equa- 
tion (A35) that, in following a constant curve, d^ will maintain a steady-state error 
after bringing to zero. 
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Vertical loop closures can be formulated in a way similar to that used for the 
lateral loop closures. The control error in this situation is based on either the vertical 
path or on the vertical direction of motion of the vehicle. The vertical path and vertical 
direction of motion are defined as the projection of the velocity vector on the ver- 

tical plane of symmetry of the vehicle. 

d v 

It has been shown in figure 38(b) that the transfer function H g (s) of the vertical, 

c 

augmented vehicle dynamics is effectively of a lower order than the corresponding trans- 
fer function of the lateral dynamics. This softens the demand on the amount of rate 
information necessary to stabilize the system. Thus, a control error based on the ver- 
tical direction might yield a satisfactory vertical loop closure, since the single zero 
might be sufficient. Furthermore, the vertical trajectory profile does not have sections 
with constant curvature, so that the problem of a steady-state error in d v does not 
occur. 
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TUNNEL DISPLAY AND GENERATION 
Tunnel Display 

The basic principle, underlying the tunnel display, is shown in figure 41. The tra- 
jectory to be followed is presented as a perspective tunnel image with a constant square 
cross section. The perspective appearance of the tunnel changes with the vehicle posi- 
tion. Positional cues are derived from this perspective appearance, while rate cues are 
derived from the rate of change in the appearance. In figure 41, the line AA represents 
the horizon and point C the center of the image, which represents the vehicle axis. 

Point C is not explicitly displayed, but the pilot can estimate its location by using the 
monitor frame as a reference. The trajectory is curved to the left, as shown in the 
corresponding plan view (fig. 42). The bank angle <p is determined from the inclination 
of the horizon with respect to the reference frame. Point F indicates the tunnel direction 
and is not explicitly displayed. However, it is assumed that the location of point F can 
be estimated by the pilot, by extending the four tunnel corner lines nearest to the vehicle, 
to a focal point at infinity. Point S is the center of the nearest square and is to be esti- 
mated as well. The centerline FS constitutes the tunnel direction at the approximate 
position of the vehicle. When defining the Fy-axis to be parallel and the F z -axis to 
be perpendicular to the horizon, the lateral and vertical deviations of the vehicle from the 
trajectory can be estimated from the distance between F and S in Fy- and F z -directions, 
respectively. 

Attitude information can be estimated from the deviation of focal point F from the 
center of the image C. Thus, the lateral intercept angle d^ of the vehicle axis with the 
tunnel direction is estimated from the distance between F and C in Fy-direction, and the 
vertical intercept angle dg of the vehicle axis with the tunnel direction is estimated 
from the distance between F and C in F z -direction. The downslope of the trajectory 0^. 
is determined by the distance point F is under the horizon AA. 

Rate cues originate from the rate of change of the above mentioned quantities. A 
relative impression of forward motion is derived from the "passing by" of tunnel seg- 
ments. Note that this velocity impression can be augmented by presenting the corner 
lines of the tunnel as separate line segments, rather than solid lines. 

A preview of the trajectory to be followed is obtained by considering a centerpoint T 
on the trajectory, a distance D ahead of the vehicle. Point R is a point on FS and is 
also a distance D ahead. Unlike R which is to be estimated, the tunnel section at T is 
explicitly shown. The distance between T and R in Fy-direction provides an indication 
of the trajectory curvature, while the distance between T and R in F z -direction indicates 
changes in the slope. 
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Predictor Display 

The lateral as well as the vertical vehicle paths are assumed to be circular and 
tangential to the velocity vector V. The lateral and vertical curvatures are given by 


1 a £ (t) 


(Bl) 


1 a v^ 
Rv v 2 


(B2) 


where a f and are the lateral (in y b body-axis direction) and vertical (in z b 
body-axis direction) accelerations with respect to the inertial reference system. The pre- 
dictor symbol P indicates a point on the predicted vehicle path, a distance D ahead of 
the vehicle. When defining C y and C z as the lateral and vertical screen axes (see 
fig. 41), the deviation of P from the center of the monitor in the Cy-direction is given by 


X 


£ “ 


v_ 

V 


+ 


1 + 


D / R t 

- ( D / R «) 2 


and in the C z direction by 


(B3) 


A v = 


w 

V 


+ 


P/Ky 

i + 0 - (d/Rv) 2 


(B4) 


where v and w are the components of V in y b - and z b -directions, respectively, and 
V is the ground speed. The values of A f and X y are in fractional screen coordinates 
(screen limits are ±1), for a ±45° horizontal-and-vertieal field of view. 

A valid control strategy for the tunnel plus predictor /director display is to mini- 
mize the error between the predictor P and tunnel section T. The linearized control 
error is then given by equation (A32) of appendix A. A prediction time Tp is chosen 
which yields the best compromise between positional accuracy and system damping for 
the given vehicle dynamics and airspeed. 

Since the vertical dynamics are of a higher bandwidth than the lateral dynamics, a 
prediction time Tp which yields the best results for the lateral control will probably be 
too large for the vertical control. An excessively large prediction time will result in a 
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sluggish response in the vertical control. Furthermore, rapid vertical predictor motions 
will occur, since the vertical disturbances are less "filtered" by the vehicle dynamics. 
This problem is solved by choosing a value of T which is the most appropriate for the 
lateral control on one hand and by reducing the vertical acceleration in computing the 
vertical curvature on the other hand. Thus, the vertical curvature is computed from 


_L = k — 

Rv P v V 2 


(B5) 


where kp y is a gain and 0 = kp V = 1. It is clear that for kp V = 0 the predictor sym- 
bol deteriorates into a velocity vector symbol. 

Where the lateral acceleration a^ is not measured or available, the curvature of 
the lateral path can be approximated from the bank angle $ according to 



Jt 

V 2 


(B6) 


Since equation (B6) only yields the correct curvature in the steady state when following 
constant curves, system performance is expected to be less good than with the predictor 
based on the lateral acceleration. However, the steady-state error in d^ in following 
constant curves with this simplified predictor will be zero as well. 


Variables Necessary To Drive the Display 

The variables necessary to display the tunnel image are the coordinates of the 
vehicle position in the inertial reference system, x i? y i? and and the attitude angles 
of the vehicle, 0/, 6, and 0. The variables x i? y^, and can only be obtained by, 
or with the help of, ground-based equipment, such as a microwave landing system. The 
accuracy with which the desired trajectory can be followed directly depends on the accu- 
racy with which the vehicle position can be estimated. The attitude angles 0 , and 0 

can be obtained with standard vertical and heading gyros. 

The variables necessary to display the predictor symbol are the velocities v, w, 
and V, as well as the accelerations a^ and a^. The accelerations a^ and can 
be measured directly with onboard accelerometers. The velocities v, w, and V can 
be estimated from the sideslip, angle of attack, and airspeed, measured at the vehicle. 
However, these variables are calculated with respect to the air mass rather than to an 
inertial reference system so that the appropriate corrections for the steady-state wind 
components have to be made. 
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An alternative method of obtaining v, w, and V is by a rotational transformation 
of Xp and z i through \p, 9, and <f). The velocities x^, y^ and can be 
derived by differentiation of x^ y^, and or derived from estimates using position 
and accelerometer data. These computations can be carried out with onboard equipment. 
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TABLE 1.- STABILITY DERIVATIVES FOR LEVEL FLIGHT 
OF CH-47C HELICOPTER 


Lateral dynamics 

Longitudinal dynamics 

Parameter 

Value 

Parameter 

Value 

Y v /m 

-0.0779 sec"* 

X u /m 

-0.0219 sec - l 

Yp/m 

-0.6246 m * sec“l 

^w/ m 

0.0358 sec"* 

Y r /m 

-0.0679 m • sec - l 

V” 

0.7146 m • sec - * 

_V m 

0.1343 m • sec - 2 cm _ l 

xh 

0.0159 m • sec“2 cm"* 

s 1 
X° i 

L 

-0.0068 m • sec - 2 cm - ^ 

xh 

0.0461 m * sec“2 cm“^ 

L v/'xx 

-0.01804 msec"l 

z u / m 

-0.0589 sec - l 

S/'xx 

-0.8210 sec"l 

Z w /m 

-0.5727 sec-1 

^ r/^xx 

-0.0675 sec"l 

Z q/ m 

-0.3357 m • sec - ^ 

hj 1 ** 

0.1597 sec -2 cm -l 

xh 

0.0562 m * sec"^ cm - * 

i 

b 

i— i 

o 

-0.0533 sec -2 cm"l 

V m 

-1.1518 m • sec~2 cm"^ 

^v/^zz 

-0.0007 (m/sec~l)“l 

Mu/^y 

-0.0138 (m/sec-l) - l 

Np/lzz 

-0.0169 sec - l 

Mw/^y 

0.0571 msec“l 

^r/*zz 

-0.0391 sec'l 

M _q/Vy 

-1.6959 sec~l 

V- 

0.0034 sec -2 cm - * 

M 5 t y/ I yy 

0.1568 sec"^ cm - 1 


0.0776 sec -2 cm"' 

M^/iyy 

0.0622 sec”2 C m"^ 

m 

14 968.6 kg 

*yy 

273 536.0 kg • m^ 

J xx 

50 386.3 kg • m 2 

U N 

33.53 m • sec”* 

*zz 

257 685.0 kg • m 2 

W N 

1.72 m • sec”l 

J xz 

19 838.3 kg- m 2 

0 N 

0.05122 rad 
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TABLE 2.- BACKGROUND OF TEST SUBJECTS 


Subject 

Occupation 

Age 

Total flight 
time, hr 

Percent of total 
flight time in 
helicopters 

Aircraft 

Previous 

simulator 

experience 

1 

Physical 

scientist 

43 

360 

15 

Cessna 150 
and 172 

Piper Cherokee 
Bell 47 

Some 


Aeronautical 

engineer 

30 

0 



Extensive 

Z 



3 

Research 

pilot 

36 

2200 

30 

All types 

Extensive 

4 

Research 

pilot 

43 

5200 

20 

All types 

Extensive 


38 



TABLE 3. - RESULTS OF TRAJECTORY FOLLOWING EXPERIMENT; SUBJECT 1 



Basic 
tunnel 
(no roll) 

EADI/map 

Tunnel plus 
vehicle axis 
cross/director 
(no roll), 

D = 305 m 

Vehicle axis 
cross/director 
(no roll), 

D = 305 m 
(1000 ft) 

Tunnel plus 
predictor /director, 

kpy = 0.2 

D = 229 m 
(750 ft) 

Tunnel plus 
simplified 
predictor 
(no roll), 

D = 229 m 

Simplified 
predictor /director 
(no roll), 

D = 229 m 
(750 ft) 




(1000 ft) 

No roll 

Roll 

(750 ft) 

No. of runs 

7 

10 

12 

10 

6 

6 

15 

3 

cov(d £ ), m2 

578.6 

3663.8 

745.5 

1840.9 

100.0 

114.8 

173.3 

189.8 

±320.5 

±2687.5 

±494.4 

±272.6 

±12.3 

±23.7 

+51.0 

±19.6 

cov^), m 2 

205.8 

1089.8 

188.5 

1775.9 

28.2 

37.8 

22.3 

21.4 

±141.7 

±746.7 

±145.6 

±605.2 

±13.6 

±10.4 

±6.8 

±3.3 

cov/cA m 2 

NR a 

NR 

4689.3 

1621.5 

146.3 

209.8 

227.0 

158.1 

V £ / 



±3092.5 

±650.9 

±27.6 

±67.3 

±102.0 

±22.6 

covfe v N ) , m 2 

NR 

NR 

412.6 

577.7 

100.0 

137.7 

105.3 

81.3 

\ V 



±97.1 

±376.3 

±29.5 

±59.2 

±42.4 

±2.2 

cov(e £ ), 

(m/sec)^ 

NR 

NR 

158.8 

80.0 

251.4 

223.3 

88.1 

60.7 



±43.9 

±19.9 

±17.1 

±20.6 

±22.9 

±6.3 

cov (e v ), 

(m/secr 

NR 

NR 

43.8 

33.3 

312.2 

117.0 

90.4 

49.9 



±13.7 

±9.0 

±13.7 

±18.5 

±39.8 

±5.6 

cov(d^), 

10’ 3 rad 2 

11.7 

51.4 

13.6 

13.9 

4.33 

4.52 

7.18 

5.17 

±3.1 

±21.7 

±3.9 

±3.9 

±.55 

±1.23 

±2.50 

±.99 

cov(p), 

8.12 

7.71 

8.69 

5.47 

2.51 

2.27 

1.63 

1.14 

(rad/sec) 2 

±2.78 

±3.89 

±3.99 

±1.32 

±.33 

±.31 

±.36 

±.07 

cov(6 g ), 

A 

2.66 

2.76 

3.03 

1.99 

1.61 

1.30 

1.16 

0.96 

±.70 

±1.29 

±1.16 

±.38 

±.18 

±.23 

±.17 

±.08 

cm z 









(roll) 









c°v(6 c ), 

10.53 

6.23 

11.92 

7.30 

6.67 

7.60 

5.99 

4.53 

cm2 

(collective) 

±4.10 

±2.59 

3.01 

±3.16 

±2.32 

±4.00 

±1.90 

±.09 


a Not relevant for this display configuration. 
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TABLE 4.- RESULTS OF TRAJECTORY FOLLOWING EXPERIMENT; SUBJECT 2 



Basic 
tunnel 
(no roll) 

EADl/map 

Tunnel plus 
vehicle axis 
cross/director 
(no roll), 

Vehicle axis 
cross/director 
(no roll), 

Tunnel plus 
predictor /director, 
D = 229 m 
(750 ft) 

Tunnel plus 
simplified 
predictor 
(no roll), 

D = 229 m 
(750 ft) 

Simplified 
predictor/director 
(no roll), 

D = 229 m 
(750 ft) 



D = 305 m 
(1000 ft) 

L/ — ovJD 111 

(1000 ft) 

No roll 

Roll 

No. of runs 

6 


12 

8 

5 

7 

4 

2 

cov ( d t)’ m2 

270.2 

1969.4 

326.4 

2485.3 

90.0 

83.6 

214.6 

155.8 

±147.0 

±1401.5 

±111.0 

±293.1 

±14.2 

±11.0 


±3.4 

cov(^d v ), m 2 

112.4 

797.6 

135.6 

1255.1 

22.9 

16.2 


25.2 

±22.7 

±538.3 

±43.0 

±311.8 

±4.1 

±3.7 


±.8 

cov (fi}’ m2 

NR a 

NR 

4947.1 

1808.0 

177.9 

202.0 

274.1 

58.0 



±975.4 

±539.3 

±21.3 

±37.6 

±57.3 

±.6 

C0V ( £ v)’ m2 

NR 

NR 

324.3 

444.0 

59.4 

mmm 

■m 

52.4 


" 

±87.2 

±243.3 

±8.2 

mSSmk 

BHkilb 

±1.2 

cov(l £ ), 

(m/sec) 2 

NR 

NR 

134.4 

148.4 

241.5 

224.5 

114.9 

105.6 



±48.7 

±66.7 

±29.4 

±20.6 

±12.1 

±8.4 

cov (*v)> 

(m/sec) 2 

NR 

NR 

22.3 

65.8 

151.1 

134.2 

59.1 

62.4 



±5.0 

±18.3 

±13.2 

±20.7 

±7.3 

±14.9 

c °v(dJ, 

10 _ 3rad2 

8.40 

38.48 

8.03 

13.55 

4.69 

4.62 

8.38 

5.13 

±1.90 

±13.25 

±2.72 

±3.36 

±.65 

±.48 

±.71 

±.21 

cov(p), 

9.30 

16.07 

16.60 

12.58 

3.62 

2.38 

2.12 

2.07 

10-3 

(rad/sec) 2 

±1.50 

±4.24 

±7.19 

±5.58 

±.41 

±.36 

±.12 

±.12 

C ° V (5 S ), 

cm2 

3.42 

5.07 

5.89 

4.12 

2.08 

1.51 

1.40 

1.63 

±.62 

±1.20 

±2.08 

±1.54 

±.23 

±.26 

±.06 

±.16 

(roll) 









c°v(6 c ), 

14.67 

8.68 

15.84 

48.64 

8.65 

8.57 

5.07 

4.92 

cm2 

(collective) 

±3.25 

±3.03 

±3.33 

±17.64 

±2.78 

±3.25 

±.41 

±.21 


a Not relevant for this display configuration. 










TABLE 5.- RESULTS OF TRAJECTORY FOLLOWING EXPERIMENT; SUBJECTS 3 AND 4 



Subject 3 


Subject 4 



Tunnel plus 
predictor /director 

Basic 

EADI/map 

Tunnel plus 
predictor /director 

Basic 



(no roll), kpy = 0.2 
D = 229 m 
(750 ft) 

tunnel 
(no roll) 

(no roll), kp V =0.2 

D = 229 m 
(750 ft) 

tunnel 
(no roll) 

EADI/map 

No. of runs 

7 

6 

7 

7 

8 

8 

cov ( d f)’ m2 

124.7 

353.2 

3850.4 

138.9 

650.3 

1652.1 

±35.2 

±286.7 

±3965.5 

±29.1 

±516.1 

±840.3 

cov^d v ), m 2 

33.7 

142.7 

313.7 

41.2 

328.8 

923.4 

±16.8 

±50.6 

±99.1 

±17.5 

±213.8 

±494.5 

cov(e £ ), m 2 

216.2 

±59.0 

NR a 

NR 

209.7 

±42.1 

NR 

NR 

cov ( e v)> m2 

131.6 

±44.3 

NR 

NR 

122.1 

±33.8 

NR 

NR 

cov(e £ ), 

(m/sec) z 

230.0 

±58.9 

NR 

NR 

239.7 
• 22.0 

NR 

NR 

c°v(e v ), 

(m/sec)^ 

115.7 

±21.0 

NR 

NR 

114.7 

±13.5 

NR 

NR 

cov ( d ,/,)> 

10 _ 3rad 2 

5.07 

12.25 

42.42 

3.66 

17.98 

37.90 

±.29 

±6.17 

±16.67 

±.96 

±7.69 

±19.32 

cov(p), 

1.63 

4.30 

3.84 

1.99 

9.82 

5.08 

10-3 

(rad /sec) 2 

±.20 

±.60 

±.52 

±.13 

±2.02 

±.87 

cov(6 s ), 

cm 2 

(roll) 

1.14 

1.66 

1.41 

1.27 

3.32 

1.94 

±.16 

±.31 

±.22 

±.10 

±.68 

±.47 

c ° v (6 ), 1 

7.64 

10.13 

5.34 

6.66 

11.79 

9.84 

cm 2 

(collective) 

±2.33 

±1.96 

±1.47 

±2.14 

±2.22 

±2.69 


a Not relevant for this display configuration. 
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TABLE 6.- RESULTS OF TRAJECTORY FOLLOWING EXPERIMENT; EFFECT 
OF D AND COMPARISON OF TUNNEL PLUS PREDICTOR/DIRECTOR 
WITH PREDICTOR/DIRECTOR; SUBJECT 1 



Tunnel plus predictor /director 
(no roll), k pv = 1 

Pr edic tor /director 
(no roll), k pv = 1 


D = 91.4 m 
(300 ft) 

D = 137 m 
(450 ft) 

D = 229 m 
(750 ft) 

D = 305 m 
(1000 ft) 

D = 91.4 m 
(300 ft) 

D = 229 m 
(750 ft) 

D = 305 m 
(1000 ft) 

No. of runs, n 

4 

5 

5 

8 

5 

5 

4 

cov^d^, m2 

*49. 2 

68.2 

107.3 

173.6 

67.0 

138.1 

276.3 

b ±5.7 

±3.7 

±26.7 

±39.6 

±4.1 

±48.5 

±41.3 

cov(d v ), m2 

11.9 

17.4 

59.2 

46.4 

12.9 

53.9 

71.0 

±1.5 

±4.2 

±40.8 

±17.7 

±2.4 

±10.5 

±10.7 

cov^, m 2 

37.4 

60.7 

260.2 

937.5 

62.1 

229.1 

515.4 

±6.9 

±4.7 

±47.7 

±102.2 

±4.8 

±54.7 

±91.7 

cov ( e v)’ m2 

23.2 

41.0 

250.3 

1 289.0 

22.0 

302.2 

1 036.2 

±2.2 

±11.2 

±88.5 

±566.2 

±2.7 

±34.0 

±129.4 

cov (^)’ 

(m/sec) 2 

13.4 

64.5 

380.7 

1 080.0 

18.2 

437.5 

694.8 

±.8 

±19.6 

±87.4 

±442.8 

±3.7 

±146.4 

±236.5 

c °v(e v ), I 

(m/sec) 2 

23.8 

160.0 

2121.7 

14 963.7 

91.6 

4897.7 

12 159.1 

±30.0 

±300.0 

±1719.3 

±4 250.6 

±8.0 

±3190.0 

1 654.0 

cov(d^), 

10'3rad2 

3.69 

4.40 

4.67 

6.98 

5.04 

4.86 

5.74 

±.46 

±.29 

±1.53 

±1.11 

±.46 

±1.85 

±1.05 

cov(p), 

3.58 

2.44 

2.24 

1.48 

4.39 

1.57 

0.89 

10-3 

(rad/sec) 2 

±.37 

±.24 

±.46 

±.37 

±.44 

±.19 

±.16 

covf 6 A 

1.68 

1.37 

1.36 

1.03 

1.79 

1.09 

0.71 

\sj> 

cm 2 

(roll) 

±.14 

±.07 

±.24 

±.20 

±.14 

±.08 

±.16 

cov(6 c ), 

5.50 

4.02 

5.59 

4.92 

5.32 

4.37 

5.54 

V W 
cm 2 

(collective) 

±.27 

±.68 

±1.07 

±1.00 

±.56 

±.39 

±1.86 


a Average of n runs. 

^Standard deviation of n runs. 
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TABLE 7.- RESULTS OF TRAJECTORY FOLLOWING EXPERIMENT; EFFECT 
OF D AND COMPARISON OF TUNNEL PLUS PREDICTOR/DIRECTOR 
WITH PREDICTOR/DIRECTOR; SUBJECT 2 



Tunnel plus predictor/director 
(no roll), kp V = 1 

Predictor /director 
(no roll), k pv = 1 


D = 91.4 m 

D = 137 m 

D = 229 m 

D = 305 m 

D = 91.4 m 

D = 229 m 

D = 305 m 


(300 ft) 

(450 ft) 

(750 ft) 

(1000 ft) 

(300 ft) 

(750 ft) 

(1000 ft) 

No. of runs, n 

4 

4 

8 

4 

4 

4 

4 

cov(d £ ), m 2 

47.4 

47.0 

102.1 

111.8 

47.8 

110.6 

227.1 

±6.2 

±3.6 

±9.9 

±11.7 

±4.4 

±14.9 

±50.3 

cov^d v ) , 

9.0 

13.5 

26.7 

33.9 

12.2 

74.9 

61.2 

±1.3 

±1.8 

±5.1 

±2.9 

±2.8 

±15.5 

±15.9 

cov(e £ ), m 2 

36.3 

42.5 

253.5 

1 272.8 

41.4 

125.8 

235.3 

±5.9 

±3.5 

±44.2 

±235.9 

±5.5 

±15.9 

±48.4 

cov(e v ), m 2 

19.6 

56.0 

342.9 

1 192.4 

22.9 

420.4 

812.2 

±1.6 

±5.2 

±47.4 

±275.4 

±5.1 

±38.5 

±11.4 

cov (€ £ ) , 
(m/sec) 2 

12.3 

76.5 

279.5 

1 144.4 

11.9 

204.3 

648.6 

±2.0 

±8.2 

±67.0 

±154.2 

±1.0 

±26.1 

±62.0 

cov(e v ), 

(m/sec) 2 

202.5 

862.9 

7570.3 

21 140.0 

420.3 

11 280.0 

20 498.0 

±95.7 

±184.1 

±3586.1 

±5 593.0 

±513.7 

±792.0 

±5 820.0 

covfdj,), 

4.18 

3.34 

4.48 

5.78 

4.38 

3.66 

3.91 

10 _ 3rad2 

±.21 

±.23 

±.74 

±.86 

±.43 

±.31 

±.40 

cov(p). 

4.14 

3.19 

2.29 

4.38 

5.13 

1.84 

1.93 

JO” 3 

±.46 

±.18 

±.30 

±.82 

±.65 

±.21 

±.26 

(rad/sec) 2 








cov ( 6 s)' 

1.89 

1.73 

1.32 

2.12 

2.06 

1.16 

1.20 

9 

±.11 

±.08 

±.14 

±.44 

±.21 

±.11 

±.22 

cm* 

(roll) 






cov(a c ), 

o 

5.85 

5.11 

4.75 

6.78 

5.81 

4.62 

5.95 

±.24 

±.43 

±.16 

±1.32 

±.33 

±.13 

±1.86 

cm z 

(collective) 
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TABLE 8.- RESULTS OF TRAJECTORY ENTRY EXPERIMENT; SUBJECT 1 



Tunnel plus predictor /director 
(no roll), k pv = 0.2 

Predictor /director 
(no roll), k pv = 0.2 

Tunnel plus 
predictor/director 
(roll), k pv = 0.2 
D = 137 m (450 ft) 

Basic 

tunnel 

(roll) 

EADI/map 

D = 91.4 m 
(300 ft) 

D = 137 m 
(450 ft) 

D = 229 m 
(750 ft) 

D = 305 m 
(1000 ft) 

D = 229 m 
(750 ft) 

D = 305 m 
(1000 ft) 

No. of series 

4 

4 

4 

4 

4 

4 

4 

4 

5 

of 6 entries 










sc(d f V m 

10.8 

11.4 

11.7 

14.8 

13.3 

20.8 

12.1 

19.4 

34.4 


±.6 

±.9 

±.7 

±.8 

±.5 

±2.5 

±1.0 

±2.1 

±3.0 

SC (dy) > m 

3.0 

3.6 

4.8 

4.7 

5.7 

7.2 

3.3 

6.9 

19.1 

±.2 

±.2 

±.3 

±-5 

±.2 

±.5 

±.2 

±.4 

±1.6 

T„, sec 

16.3 

16.8 

16.2 

22.1 

23.4 

36.0 

17.1 

28.3 

46.5 


±.6 

±1.0 

±.4 

±4.0 

±1.4 

±1.3 

±1.4 

±5.6 

±2.9 

T , sec 

10.2 

10.5 

12.9 

14.6 

18.1 

24.2 

10.6 

31.3 

44.6 


±.5 

±.6 

±2.3 

±2.4 

±1.7 

±3.7 

±.2 

±2.1 

±2.3 

cov(e 

64.0 

48.0 

21.6 

17.0 

20.3 

11.1 

47.9 

NR a 

NR 

V xj 

±1.3 

±1.1 

±.5 

±2.0 

±.5 

±3.5 

±1.9 



10 2 m 2 










cov(e v ) 

1.20 

0.91 

0.80 

1.30 

1.17 

0.91 

0.94 

NR 

NR 

\ v / 

O n 

±.10 

±.04 

±.06 

±.20 

±.14 

±.19 

±.06 



10 2 m 2 










C0V(€.), 

1.03 

0.94 

1.17 

2.64 

1.11 

2.56 

0.99 

NR 

NR 

o 9 

±.03 

±.04 

±.18 

±.88 

±.07 

±.67 

±.04 



10 2 (m/secr 










cov(e v ), 

0.09 

0.22 

0.91 

1.96 

0.97 

2.97 

0.20 

NR 

NR 

10 2 (m/sec) 2 

±.01 

±.01 

±.33 

±.38 

±.63 

±.69 

±.04 



cov(p), 

3.39 

2.37 

1.96 

1.91 

1.96 

1.57 

2.72 

6.02 

9.40 

10-3 

±.15 

±.23 

±.27 

±.23 

±.15 

±.17 

±.33 

±1.55 

±2.36 

(rad/sec) 2 










c°v(6 ), 

1.26 

0.94 

0.85 

0.81 

0.86 

0.70 

1.15 

1.94 

3.15 


±.05 

±.09 

±.10 

±.09 

±.06 

±.05 

±.12 

±.48 

±.73 

cm 2 (roll) 










c°v(6 c V 

2.92 

2.72 

2.46 

2.50 

2.21 

2.08 

2.52 

5.69 

4.39 

\ / 
9 

±.10 

±.05 

±.04 

±.21 

±.09 

±.06 

±.10 

±.75 

±1.71 


cm2 

(collective) 


a Not relevant for this display configuration. 























TABLE 10.- RESULTS OF TRAJECTORY ENTRY EXPERIMENT; SUBJECT 3 



Tunnel plus predictor /director 

Predictor /director 




(no roll), k nv — 

0.2 

(no roll), 

r 

II 

O 

to 

Basic 





— 



tunnel 

EADI/map 


D = 137 m 

D = 229 m 

D = 305 m 

D = 229 m 

D = 305 m 

(roll) 



(450 ft) 

(750 ft) 

(1000 ft) 

(750 ft) 

(1000 ft) 



No. of series 

3 

3 

3 

3 

3 

3 

4 

of 6 entries 








sc(d £ ), m 

11.3 

16.9 

17.8 

22.5 

19.7 

22.7 

39.7 

±.6 

±2.2 

±1.9 

±3.0 

±1.3 

±.9 

±14.0 

sc(d v ), m 

5.27 

6.73 

8.45 

9.90 

11.61 

6.82 

22.0 

±.18 

±.50 

±.90 

±1.15 

±.81 

±.84 

±8.4 

T*, sec 

18.0 

19.4 

34.2 

28.0 

40.7 

29.0 

48.3 

SJt 

±.6 

±.6 

±5.3 

±.6 

±3.6 

±.2 

±1.4 

T , sec 

12.9 

16.7 

25.4 

28.6 

32.3 

27.4 

41.4 

bv 

±1.1 

±3.0 

±2.6 

±3.1 

±2.4 

±3.5 

±4.8 

covfe ), 

47.6 

26.7 

26.2 

30.5 

9.0 

NR a 

NR 

V * J 

±1.3 

±.8 

±8.7 

±5.6 

±.5 



10 2 m 2 








cov(eJ\, 

1.37 

1.53 

2.62 

2.61 

2.87 

NR 

NR 


±.05 

±.18 

±.67 

±.54 

±.53 



10 2 m 2 








covfeA 

1.05 

0.87 

4.28 

1.43 

1.82 

NR 

NR 

10 2 (m/sec) 2 

±.03 

±.08 

±1.80 

±.14 

±.14 



c°v(e v ), 

0.17 

0.47 

0.91 

0.15 

0.27 

NR 

NR 

±.04 

±.21 

±.50 

±.02 

±.06 



10 2 (m/sec) 2 








cov(p), 

2.31 

1.35 

3.18 

2.24 

1.78 

5.51 

4.36 

10-2 

±.14 

±.19 

±1.18 

±.23 

±.02 

±.65 

±.63 

(rad/sec)^ 








cov ( 6 s)’ 

0.92 

0.58 

1.11 

0.88 

0.70 

1.74 

1.48 


±.05 

±.07 

±.36 

±.07 

±.01 

±.18 

±.22 

cm 2 (roll) 








cov ( 6 c)> 

2.66 

2.10 

2.45 

1.72 

1.68 

5.96 

1.75 

\ w 

9 

±.14 

±.11 

±.30 

±.08 

±.02 

±.94 

±.29 

cm z 








(collective) 









a Not relevant for this display configuration. 
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TABLE 11.- RESULTS OF TRAJECTORY ENTRY EXPERIMENT; SUBJECT 4 



Tunnel plus predictor /director 
(no roll), kp V = 0.2 

Predictor /director 
(no roll), kp V = 0.2 

Basic 

tunnel 

(roll) 

EADI/map 


D = 137 m 
(450 ft) 

D = 229 m 
(750 ft) 

D = 305 m 
(1000 ft) 

D = 305 m 
(1000 ft) 

No. of series 
of 6 entries 

2 

4 

4 

3 

3 

4 

sc(d £ ), m 

17.5 

±2.0 

18.7 

±2.1 

24.5 

±2.2 

26.8 

±2.2 

25.5 

±4.5 

39.7 

±14.0 

sc (dy^ > rn 

4.26 

±.48 

5.33 

±.28 

6.11 

±.27 

8.73 

±.53 

12.05 

±1.69 

22.00 

±8.26 

T sr sec 

20.6 

±1.4 

24.4 

±2.0 

38.8 

±2.9 

36.8 

±2.2 

42.5 

±2.0 

48.3 

±1.4 

T sy , sec 

12.0 

±1.5 

15.4 

±3.3 

19.0 

±4.0 

27.6 

±1.8 

42.8 

±1.4 

41.4 

±4.8 

cov(e £ ), 
10 2 m 2 

55.4 

±.3 

25.5 

±1.6 

21.9 

±7.0 

16.3 

±3.0 

NR a 

NR 

cov ( e v) » 
10 2 m 2 

1.29 

±.08 

1.25 

±.36 

1.43 

±.32 

1.87 

±.31 

NR 

NR 

cov(e £ ), 

10 2 (m/sec) 2 

1.15 

±.02 

1.16 

±.05 

3.26 

±1.26 

3.55 

±2.15 

NR 

NR 

cov (e v ), 

10 2 (m/sec) 2 

0.19 

±.01 

0.88 

±.54 

2.26 

±1.24 

0.79 

±.37 

NR 

NR 

cov(p), 

10-3 

(rad/sec)* 4 

2.60 

±.37 

1.78 

±.15 

2.22 

±.74 

1.27 

±.10 

9.30 

±.34 

4.36 

±.63 

cov(6s), 
cm 2 (roll) 

0.95 

±.10 

0.72 

±.05 

0.83 

±.27 

0.54 

±.04 

2.95 

±.18 

1.48 

±.22 

cov(5 c ), 

cm 2 

(collective) 

2.74 

±.11 

2.74 

±.53 

2.35 

±.19 

1.97 

±.20 

8.49 

±2.48 

1.75 

±.29 


a Not relevant for this display configuration. 
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TABLE 12.- RESULTS OF MONITORING EXPERIMENT; SUBJECTS 1 AND 2 



Subject 1 


Subject 2 



Tunnel plus predictor /director 

Predictor/director 

EADl/map 

Tunnel plus predictor /director 

Predictor /director 



(no roll), k pv = 0.2 
D = 229 m (750 ft) 

(no roll), k pv = 0.2 
D = 229 m (750 ft) 

(no roll), k pv =0.2 
D = 229 m (750 ft) 

(no roll), kpy = 0.2 
D = 229 m (750 ft) 

EADI/map 

No. of series 

13 

12 

12 

19 

18 

16 

t of 4 







* T r , sec 

8.3 

9.0 

11.2 

7.5 

8.1 

8.6 

±1.8 

±1.6 

±1.5 

±1.4 

±1.6 

±1.5 

T c , sec 

6.6 

7.3 

10.7 

5.0 

9.0 

8.8 

±1.3 

±2.3 

±4.7 

±2.8 

±2.4 

±3.6 

eovfd.V 

2.04 

3.95 

9.29 

1.62 

3.81 

j 4.88 ' 

10 2 m 2 

±.58 

±2.87 

i 

±9.29 

±.26 

±1.77 

! ±5.03 

■ cov(d v ), 
10 2 m 2 

0.67 

i 0.90 

1.49 

0.60 

0.84 

0.96 

±.14 

| 

±.37 

i 

±.48 

±.12 

i 

±.26 

l 

; ±-53 | 

cov ( € f)’ 
10 2 m 2 

| 1.72 

±.59 

i 

i 

2.69 

±1.66 

NR a 

| 1.68 

| ±.63 

3.16 

±1.57 

NR 

C(WeA 

o o 

0.87 

±.21 

1.03 

±.39 

NR 

1 0.68 
±.17 

0.84 

±.23 

NR 

10* m ! 

c°v(d i// ), 

9.70 

14.47 

27.65 

6.83 

12.44 

20.07 

±3.20 

! ±9.39 

±16.32 

±1.76 

±5.65 

±9.63 

lO^rad* 

cov(p), 

2.25 

2.36 

6.75 

i 2.39 

■ 1.94 

8.00 

10' 3 

±.34 

±.62 

±3.36 

; ±.46 

±.34 

±2.59 

(rad /sec) 2 1 

c°v(6 s ), 

1.44 

1.58 

2.96 

‘ 1.46 

1.38 

3.15 

\ ° / 

±.15 

±.35 

±1.05 

±.18 

±.26 

±.75 

cm^ (roll) 

cov(6 c ), 

7.63 

6.99 

9.15 

7.38 

8.15 

7.72 

±1.95 

; ±2.07 

±2.06 

±1.74 

±2.03 

±1.97 

cm z 

, (collective) 

-L 

1 

| 

1 

j 

i 

i 



a Not relevant for this display configuration. 
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Figure 1.- Conceptual tunnel display. 
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(c) Predictor /director display. 
Figure 2.- Continued. 
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(e) Vehicle axis cross and director square display. 
Figure 2.- Continued. 
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Figure 3.- Functional diagram of experimental system, 
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Initial vertical deviation above the trajectory 
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Yi (East) , m (ft) 

Figure 5.- Plan view of desired trajectory. 



Traveled distance along trajectory, m 

Figure 6.- Descent profile of desired trajectory 
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director 
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Figure 8.- Results of trajectory following experiment; vertical deviation score, 

subject 1. 
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(a) Tunnel/predictor display, 
k pv = 0.2, 

D = 305 m (1000 ft). 


Figure 10.- Time histories of 



(b) Basic tunnel display. (c) E ADI/ map, 


response for trajectory following; subject 1 
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(a) Tunnel/predictor display, (b) Basic tunnel display. (c) EADI/map. 

kpv = 0.2, 

D = 305 m (1000 ft). 


Figure 11.- Time histories of vertical response for trajectory following; subject 1. 
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Figure 12. - Results of trajectory following experiment; effect of D on lateral 

deviation score, subject 1. 
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Predictor distance, D, m(ft) 

Figure 14. - Results of trajectory following experiment; effect of D on 

roll activity, subject 1. 
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(a) Tunnel plus predictor /director; (b) Tunnel plus predictor/director 
kpv = ■*’ kp v = 

D = 91.4 m (300 ft). D = 229 m (750 ft). 

Figure 15.- Time histories of lateral response for trajectory following; 

effect of D, subject 1. 


m 






kpy - 1> kp V - 1, 

D = 91.4 m (300 ft). D = 229 m (750 ft). 


Figure 16.- Time histories of vertical response for trajectory following; 

effect of D, subject 1. 



(a) Tunnel plus predictor /director; (b) Predictor/director; 

D = 305 m (1000 ft). D = 305 m (1000 ft). 


Figure 17.- Time histories of lateral response for trajectory following; comparison 
of tunnel plus predictor /director with predictor /director, subject 1. 
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(a) Tunnel plus predictor/director; (b) Predictor/director; 

D = 305 m (1000 ft). D = 305 m (1000 ft). 

18.- Time histories of vertical response for trajectory following; comparison 
of tunnel plus predictor /director with predictor/director, subject 1. 
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Figure 19.- Results of trajectory entry experiment; 
lateral deviation score, subject 1. 
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Figure 20.- Results of trajectory entry experiment; 
vertical deviation score, subject 1. 
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Figure 21.- Results of trajectory entry experiment; 
lateral settling time, subject 1. 
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Figure 22.- Results of trajectory entry experiment; 
vertical settling time, subject 1. 
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(a) Tunnel display. (b) EADl/map display. 


Figure 24.- Time histories of lateral response for entry maneuvers; comparison 
of basic tunnel display with EADl/map, subject 1, initial location 3. 
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(a) Tunnel display. 


(b) EADI/map display. 


Figure 25.- Time histories of vertical response for entry maneuvers; comparison 
of basic tunnel display with EADI/map, subject 1, initial location 3. 
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Figure 26.- Results of trajectory entry experiment; effect of D on lateral 

deviation score, subject 1. 
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Predictor distance, D, ra(ft) 

Figure 27.- Results of trajectory entry experiment; effect of D on vertical 

deviation score, subject 1. 
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Figure 28.- Results of trajectory entry experiment; effect of D on lateral 

settling time, subject 1. 
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Figure 30.- Results of trajectory entry experiment; effect of D on roll activity, 

subject 1. 
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(a) Tunnel plus (b) Predictor /director; 

predictor /director; D = 305 m (1000 ft). 

D = 305 m (1000 ft). 
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Figure 33.- Time histories of lateral response for entry maneuvers; comparison 
of tunnel plus predictor /director with predictor/director; initial location 4. 
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(a) Tunnel plus (b) Predictor /director; 

predictor /director; D = 305 m (1000 ft). 

D = 305 m (1000 ft). 


Figure 34.- Time histories of vertical response for entry maneuvers; comparison of 
tunnel plus predictor /director with predictor /director; initial location 4. 
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(a) SAS disengaged. 



(b) SAS engaged. 


Figure 37.- Complex-plane plot of pole-zero location of (s) 


j = imaginary axis, w = frequency (rad/sec), and <? = t>’. r 









(b) Based on vehicle velocity vector. 
Figure 40.- Block diagram of vehicular control. 
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